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ABSTRACT 

A forecast model proposed by 'Arnason as being capable of producing 
realistic forecasts of displacement of systems at and below 500 mb is 
investigated. The prediction equation, 

IS -V7--9| 4F (2 WWvP) as 
Jt moe \Ja 

is applied to test cases at 1000, 850, 700 and 500 mb with various values 
of the parameter k;}. The effect of the value of the parameter on the 
amount of divergence present in the model is discussed. Root-mean-sgquare 
errors for 2l-hour and )\8-hour forecasts with various values of k, at 
each level are presented. Optimum values of ky for the cases tested are 
selected as 3.5, 2.75 and 0.9 for the levels 850, 700 and 500 mb respec- 
tively. It is found that with proper values for ky the model will provide 
good forecasts of system movement at each level tested, but as expected, 
changes in the intensity of systems are not correctly forecast. In all 
cases values of root-mean-square error of forecasts are less than those 
obtained by applying the "Helmholtz" barotropic model. A geostrophic 
version of the model is also tested and is found to produce results simi- 
lar to those obtained for the stream-function version except that the 
geostrophic "blow-up" of low-latitude high pressure cells is observed. 

The writer wishes to express his appreciation for the guidance and 
encouragement given him by Professor George J. Haltiner of the U. S. 
Naval Postgraduate School in the course of this investigation. In addi- 
tion, appreciation is expressed to Mr. Geirmundur ‘'Arnason, on whose work 
this paper is based, and to personnel of the Fleet Numerical Weather 


Facility for many instances of assistance. 
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ieeewnereduction 

Since the initation of operational numerical weather forecasts by 
the Joint Numerical Weather Prediction Unit in 1955, the only forecast 
model which has stood the test of daily operational use has been the 
barotropic model first proposed by Rossby in 1939. The barotropic model 
remained in its original form until 1958 when Wolff [1] discovered large 
systematic errors due to improper forecast of very long waves. After 
these errors had been empirically corrected by Wolff, it was discovered 
by Cressman [2] that these long-wave errors could be reduced by the in-= 
troduction of divergence into the barotropic model. This resulted in 
the divergent or "Helmholtz" version of the barotropic model which is in 
general use at this time. | 

With the improved barotropic model, high quality forecasts have been 
available for the 500-mb level. This barotropic model, however, is by 
nature restricted in its use to the 500-mb level and, in addition, is ca= 
pable mainly of forecasting the movement of systems, and not in general 
changes in intensity. Occasionally attempts have been made to replace 
the barotropic model in operational use with some type of baroclinic 
model which would obviate these shortcomings. These attempts have so far 
met with Satine largely due to the inability to control the forecast of 
cyclogenesis. 

‘Arnason has suggested that, in view of the great utility of the 
500-mb forecasts produced by the barotropic model, it might be worth- 
while temporarily to set aside the problem of predicting cyclogenesis 
and to try to develop a model which would possess the same skill of fore- 
casting movement at other levels which the barotropic has at 500 mb. To 
this end, 'Arnason has proposed a model [3] which maintains many of the 
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characteristics of the barotropic, but which is more adaptable to other 
levels. 

This model is based on a single-layer compressible fluid with a 
stratified density distribution, bounded by two constant-pressure sur- 
faces. The theoretical development results in an expression for the 


horizontal wind divergence 


dus 1 (ae re) (1) 
3° H ol t 


Introducing this expression into the vorticity equation, 


IS t+ WV V7-adw V Os (2) 
by 


gives the prediction equation for the "stratified" model 


YS oe WVI9-4Y_ [IZ eWVP\z 0. (3) 
~~ es an 2 


In order to apply this equation to forecasting, 'Arnason has suggested, 
as one alternative, the use of a stream function for computation of wind 
and vorticity and the elimination of d through a geostrophic assump-= 


tion. This gives 


VCH \VIV eWVY-E% VV V=O (4) 
( ot) (FR 


which is of the same form as the prediction equation for the divergent 
or "Helmholtz" barotropic model, but with an additional divergence term 
‘namtetersoare WY) and with the substitution of the stability-dependent term 
Ee or -<. 

Although the forms of the prediction equations are similar, there 
are several basic differences between the "stratified" and other barotrop- 
ic models. By the method of perturbation analysis 'Arnason has obtained 
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the meteorologically significant root of the frequency equation for the 


stratified model 


(5) 
C Guise PAS 


Although this equation is again similar to the corresponding equation 
for the Helmholtz barotropic model, an important difference arises in 
the fact that the basic current U does not have a multiplying factor. 
By empirically varying the value of g'H, the magnitude of the @ term 
may be changed without reducing the effect of the basic current, thus 
allowing control over the speed of movement of systems. As can be seen 
from equation (1) this is equivalent to a variation in the amount of 
divergence present in the model. This provides a method of adjusting 
the model to yield proper system movement at a number of levels. 

The stratified model as proposed by 'Arnason seems to provide a 
promising approach to numerical forecasting at and below the 500=mb 
level. It is the purpose of this paper to apply the model to a number 


of test cases and determine its effectiveness. 





2. Application of the Model to Forecasting 

In the application of the prediction equation (3) to actual fore- 
casts there are several choices to be made. The parameter V has not 
been explicitly defined, but may be taken to be any kind of space- 
averaged or zonal wind. Also the value of g'H is left to be determined 
empirically. 

By examination of the prediction equation (3), it can be seen 
that the additional "divergence" term of the stratified model which in= 
volves V essentially represents the horizontal advection of height 
by the mean wind. This immediately suggests the notion of a steering 
current, which may be thought of as a wind at some level or as a mean 
wind in some layer of the atmosphere. Also, the magnitude of the di- 
vergence term is dependent on the difference between the mean wind and 
the wind at the forecast level. Thus it appears reasonable to use as 
a mean wind something close to the 500=-mb wind, since this would mini- 
mize the divergence term at 500 mb and allow it to increase with dis- 
tance from 500 mb. This is generally consistent both with observed 
patterns of divergence in the atmosphere and with the need for increase- 
ing the value of the divergence term at lower levels to supplement the 
movement due to vorticity advection alone. 

A second parameter which must be empirically interpreted is the 
factor g'H. The value of g' is a measure of the static stability of 
the basic flow and H is the thickness of the layer considered in the 
model. Since neither is uniquely defined, the combination may be taken 
as any reasonable value which produces proper movement of the systems 
being forecast. 

In order to make a forecast, either Y or Pf must be eliminated 


hi 





from equation (3). ‘'Arnason has used a geostrophic relation to eliminate 
p , obtaining equation (li) which involves alone. To apply this equa- 
tion, height fields must first be "balanced" to obtain the stream function 
and then, at the conclusion of the forecast, inverted to obtain the fore-~ 
cast height field. Another approach would be to compute the vorticity 

and winds directly from the height values using the geostrophic rela- 
tionship, obtaining the forecast equation: 


V8) ie + SV V7 > We V2 = © (6) 
( on) 3 3 / g" 





3. Procedure 

As mentioned above, in order to apply the stratified model to actual 
forecasting, values must be chosen for the parameters WV and g'H. For 
the purposes of the tests presented here, it was decided to use Aas 
a mean wind computed from a height field obtained by first heavily 
smoothing the 500=-mb and 700=mb heights, then computing an intermediate 


field by the equation 


$= 715 +.25¢_., (7) 


200 TOO 


The wind obtained by this method is thus a space-mean wind for some 
intermediate level between 500 mb and 700 mb. For the first six one- 
hour time steps, a WV computed from the initial data was used; there- 
after a new WV was computed from forecast fields every sixth time step. 

In the choice of the parameter g'H, consideration was given to the 
similarity between the forecast equation (4) and the forecast equation 
for the "Helmholtz" barotropic 

(V KEY ) 2 tVWVD=0, (8) 
ie | et 

Since it was known that forecasts made with this model at 500 mb are of 
consistently good quality, it was decided to choose g'H such that the 
Helmholtz term for the stratified model would be the same as the Helm- 
holtz term of the barotropic. Thus in the Helmholtz term, g'H was 


chosen to be 
gl = )ihusin (9) 
ke 
where k is a "tuning constant" with a value of four, and Hp is the stand- 
ard height of the 500-mb surface. 
In order that the model be used at levels other than 500 mb, it was 
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decided to include a second "tuning constant" in the final divergence 
term. By varying this second constant the prcper amount of divergence 
may be introduced at each level to produce the best forecast. With 
these values introduced, the forecast equation for the stream function 
case becomes 


( v-4E7 Maw Vth ta Vvreo (10) 
7) / qH< 


In the same manner, starting with equation (6) a similar equation was 


obtained for the geostrophic version of the model: 


(v-3F 2) J? tEWV 7 ~ EE? V:V2Z=0 (11) 
Gis / dt 38 Gus 


The model, in the form of finite difference equations obtained from 
equations (10) and (11), was run to }\6-hour forecasts for 850 mb, 700 mb 
and 500 mb for several randomly selected days. Various values of k, were 
used at each level for both the geostrophic and stream-function cases. 

In addition, for the 6 May 1962 case runs were made on 1000-mb data in 


order to determine the feasibility of using the model at that level. 





4. Computational Methods 

The test computations of the stratified model described in this 
paper were made on a Control Data Corporation 160); computer using a 
program constructed of meteorological subroutines prepared by the Fleet 
Numerical Weather Facility, Monterey. The grid used in the numerical 
solution of the prediction equation was the standard octagonal grid used 
by both the Joint Numerical Weather Prediction Unit at Suitland, Md. and 
the Fleet Numerical Weather Facility. 

Time differencing was accomplished in one-hour time steps using the 


central difference formula 
Po. . ‘ T= tr 2 ( At), ) (12) 


where A, ils the computed one-hour change in geopotential. At the 
initial time, Y=O , and all integral multiples of twelve, the forward 
difference formula 

Seog = Paes (A, >). a) 
was used. 

In cases where a stream function wind was used, the stream values 
were obtained from initial height data through the solution of the bal- 
ance equation as described by 'Arnason fy). All initial height fields 
and fields required for forecast verification were obtained from files 


of operational analyses at the Fleet Numerical Weather Facility. 





5. Results 
Verification scores for each of the cases tested were determined 


by a FNWF program using the equations 
N 


Pillow = iS ( A-B),, (1h) 
N 


| N 2 
RMSE = = Lie oiler. (15) 
N 


V 


Here A represents the forecast value, B, the observed value and N, the 
total number of grid points, 1977. | 

As can be seen from the equations, the pillow is a measure of the 
error in the forecast of the fees height while the RMSE (root mean square 
error) is a measure of the deviation of errors around the mean error. 

Figures 1, 2 and 3 show plots of RMSE as a function of ky for the 
various cases tested. As can be seen, the value of k, necessary to pro- 
duce the best forecast (minimum RMSE) varies somewhat from case to case 
at each level. However, the curves are relatively flat in the vicinity 
of the minima so that by selecting a reasonable value of kas the qualiege 
ty of forecasts from one case to another will vary only slightly from 
the optimum. Table 1 shows mean values of RMSE for all cases obtained 
with various values of ky at each level. A plot of these values is 
shown in figure ). From this plot of mean error versus ki> values of 
k, can be selected which will give the minimum total error. Table 2 
gives a tabulation of the error which would occur in each case tested if 
k, at each level were selected as the value which produces the smallest 
mean error. 

Forecasts for 850 mb made from initial data for 20 January 1962 are 
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Table 1. Mean Verification Scores of Forecasts With Various Values of ky 


ee a ee | 


STREAM=-FUNCTION FORECASTS 


ee ee 


RMSE (ft) |} k, | RMSE (ft) }] ky | k, | RMSE (ft) 





Ky 
| 850 2 Tare 3 107 108 5 11h 
S, ; 
S| 700 | 1] 133 2| 12h 123 bh} 128 
: S| soo O 14:8 > 147 17 2 150 
| 
: ee 
|»; 850 | 2] 1h 3} 166 h| 166 5} 173 
tg 
{ 
| g 700 i 189 2 174 2 169 rn 173 
SS 500 0 219 $ 215 1 21) 2 218 
| GEOSTROPHIC FORECASTS 
: ky | RMSE (ft) || ky | RMSE (ft) || ky | RMSE (ft) || ky | RMSE (ft) 
— ° i 
mh oee 3 123 ly 125 5 13) 
= 
8 700 2 127 3 126 \ 131 
| S| 500 =) 461 1| 162 2| 169 
| 
| » | 850 3 192 \ 193 5 200 
e 
: g 700 2 203 3 195 198 
| =i lll Bao > Fil 1 269 2 277 








Table 2. Verification Scores For Each Day Using Optimum Value of kj 








850 mb 700 mb S00 mb 
RMSE (ft) RMSE (ft) RMSE (ft) 
ky =3.5 ki 40. 9 
), November 1961 ULL 116 13h 
(2h hr) 
20 January 1962 128 Ze 178 
(2h hr) 
6 May 1962 87 88 127 
(2) hr) 
h November 1961 173 185 
(48 hr) 
20 January 1962 198 185 250 
(48 hr) 
6 May 1962 ns 1h0 168 


(48 hr) 





shown in figures 8 through 13. Initial charts and charts for the veri-~ 
fying time of each forecast are shown in figures 5 through 7. Comparing 
forecasts for the various values of k, at a given level shows that the 
forecast speed of synoptic waves increases as the value of ky is increased 
(thus decreasing the value of g'H and increasing the magnitude of diver- 
gence). This is in accordance with the results of equation (5). Further, 
it can be seen that there is very little change in intensity of the vari- 
ous systems as k; is varied. This demonstrates the point that the addi- 
tional divergence term of the stratified model results in movement rather 
than development. 

As can be seen from equations (10) and (11) a choice of k=0 would 
reduce the forecast to the veer An examination of figure } re- 
veals that the optimum value of k, at 500 mb is close to zero, but that 
the stratified model does show a slight improvement over the barotropic 
with values of k, between zero and one. This is as expected, since 500 
mb is generally accepted to be at or near the level of minimum divergence, 
and thus the contribution of divergence to the change in the vorticity 
pattern at this level should be small. 

At levels other than 500 mb, where the barotropic model is less ade- 
quate, the values of k, for which the minimum RMSE is reached are greater. 
For the 700-mb level optimum values for the cases tested range from ee 
to k1}=, while at 850 mb they range from 3 to 5. These results can be 
accounted for by the vertical distribution of velocity divergence which 
generally increases with distance from the 500-mb surface. 

By use of equation (1) and the chosen values for g'H given by equa- 
tion (9), an expression for the divergence at a given time may be ob- 
tained: 
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AGN ie 0G Gr: + ie (16) 
dst \ da 


In figures 1), 15 and 16 divergence charts at 850, 700 and 500 mb as com- 
puted from equation (16) for initial data of 20 January 1962 are shown. 
These charts show clearly the vertical distribution of divergence, with 
values as large as + 5 x 107© sec” at 850 mb, + 3x 107° soe. at 700 
mb and values very nearly zero at 500 mb. 

An examination of the divergence fields relative to their correspond- 
ing height fields reveals that lines of zero divergence very closely fol- 
low the short-wave trough and ridge lines of the corresponding height 
fields. Between the sheet uae trough and ridge lines are located the 
centers of maximum convergence and divergence, with divergence ahead of 
ridges and convergence ahead of troughs, as might be expected, 

For the case of 6 May 1962, 1000-mb forecasts were made. From 
figure 3 it can be seen that the curves of RMSE versus k, for 1000 mb 
are almost exactly parallel to the curves for 850 mb with an RMSE about 
ten feet greater at 1000 mb than at 850 mb. The fact that the curves 
are nearly parallel suggests that the WEEE cIoe” error at 1000 mb is due 
to a greater change in the intensity of systems at this level than at 
850 mb. The movement of established systems which appear on both initial 
charts is nearly the same at both levels; thus k, which controls the mag~ 
nitude of the additional movement term should reasonably be chosen the 
same for the two levels. The verifying charts for the forecasts in 
question, figures 21, 22, 25 and 26, show this to be the case. Forecasts 
at both levels have the major systems placed very close to their actual 
positions at verifying time, but the verifying charts at 1000 mb differ 
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more in pattern due to the variation in intensity change of different 
systems over the forecast period. 

In addition to the forecasts made for the three days using stream 
functions obtained from the initial height fields, for two of the days 
forecasts were made directly from the height fields using equation (11). 
As can be seen from figures 27 through 30, the geostrophic forecasts 
thus obtained showed great similarity to the forecasts made from the 
stream function. However, as expected, the geostrophic version appears 
to suffer from the major error common to geostrophic forecasts = the 
"blow-up" of low-latitude high-pressure areas. As can be seen, the 6- 
hour geostrophic forecast increases the central height of the 500=mb 
high off the east coast of enetunitad States to over 19,500 feet, while ° 
the stream-function version keeps the central height to a more reasonable 
value. The treatment of the high off the west coast of the United States 


is similar. 


Lh 





6. Conclusions 
The stratified model as proposed by 'Arnason appears, on the basis 


of the test cases presented here, to produce reasonable forecasts for 
levels at and below 500 mb. The stream-function version of the strati- 
fied model produces 500-mb forecasts for 2); and 8 hours which are in 
all cases very slightly better than those produced by the barotropic 
model. At 850 mb and 700 mb the model produces forecasts of system 
movement which are very realistic. The geostrophic version of the model 
suffers from the "blow-up" of low-latitude high pressure cells, as is 
common With most geostrophic models. 

Optimum values of the multiplying factor of the divergence term 
k, at 850, 700 and 500 mb have been found to be 3.5, 2.75 and 0.9 
respectively. Although the number of cases tested was not large, the 
values determined here should be sufficiently accurate to provide a 
basis for the long-term testing necessary before any model becomes suite 


able for operational use. 
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Figure 1. Plot of k, vs. RMSZ for 850 mb 
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Figure 3. Plot of ky 
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Figure . Plot of k, vs. mean RMSE 
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Figure 5. 850-mb analysis, 20 January 1962 
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Figure 6. 850-mb analysis, 21 January 1962 








Figure 7. 850-mb analysis, 22 January 1962 
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Figure 8. 2lj-hour forecast for 850 mb from 20 January 1962, k, =0 
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Figure 9. 2-hour forecast for 850 mb from 20 January 1962, 








Figure 10. 2h-hour forecast for 850 mb from 20 January 1962, k,=5 
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Figure 11. ),8-hour forecast for 850 mb from 20 January 1962 
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Figure 12. 48-hour forecast for 850 mb from 20 January 1962, k,=3 
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13. 48-hour forecast for 850 mb from 20 January 1962, k,=5 
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Figure 15. 700-mb divergence, 20 January 1962 
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Figure 16. 500-mb divergence, 20 January 1962 
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Figure 18. 850-mb analysis, 6 May 1962 
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Figure 19. 2h-hour forecast for 1000 mb from 6 May 1962 
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Figure 20. 2h-hour forecast for 850 mb from 6 May 1962 
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1000-mb analysis, 7 May 1962 
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Figure 21. 








Figure 22 
- 850-mb analysis, 7 May 1962 


37 








y 1962 


4a 


hour forecast for 1000 mb from 6 ! 


Figure 23. l8- 





hs 


(Aeron 
Ay eS 


¥ 


Leer 





Figure 2). 8-hour forecast for 850 mb from 6 May 1962 
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Figure 26. 850-mb analysis, 8 May 1962 








Figure 27. 500-mb analysis, 4 November 1961 
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l8-hour forecast for 500 mb from November 1961 
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Figure 28. 
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h8-hour forecast for 500 mb from ); November 1961 
( geostrophic ) 


Figure 29. 








Figure 30. 500-mb analysis, 6 November 1961 
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